Interaction regimes for oppositely charged plates with multivalent counterions 
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Within a mean field treatment of the interaction between two oppositely charged plates in a salt 
free solution, the distance at which a transition from an attractive to a repulsive regime appears can 
be computed analytically. The mean field description however breaks down under strong coulombic 
couplings, that can be achieved at room temperature with multivalent counter-ions and highly 
charged surfaces. Making use of the contact theorem and simple physical arguments, we propose 
explicit expressions for the equation of state in several situations at short distances. The possibility of 
Bjerrum pair formation is addressed and is shown to have profound consequences on the interactions. 
To complete the picture, we finally consider the large distance limit, from which schematic phase 
diagrams discriminating attractive from repulsive regions can be proposed. 

PACS numbers: 



I. INTRODUCTION 

T ^> 1 without pair contribution Although it has been 
less studied than its like-charge counterpart [H-Q , the be- 
haviour of two interacting oppositely charged mesoscopic 
bodies in solution is of importance in various contexts, 
including colloid physics Wa|, biochemistry related ex- 
periment interpretations [oL Il0j|. drug design and 
structural biology [l2| . The simple system of two charged 
plates with opposite uniform surface charges represents 
a model of choice that enables one to get analytical re- 
sults in some limits, and furthermore, that provides a 
starting point to estimate the interaction energy between 
two colloids of various geometry [ID, [lij]. It has been 
shown that opposite charge re pul sion could occur within 
a mean field (MF) treatment [l5l [l6j . The physical ori- 
gin of such a repulsion has been identified to be twofold: 
a Born repulsion due to short range polarization effects 
when the solvent has a dielectric constant that signifi- 
cantly exceeds those of the macromolecules [3, [l7[ , and 
an osmotic repulsion resulting from the trapping of the 
counter-ions, that ensure electroneutrality between the 
unequally charged plates [H, Il9j . Within a mean-field 
approach for q : q symmetric solutions (with q = 1), it 
has been also emphasized recently that the osmotic re- 
pulsion may explain how proteins' shape determines their 
interaction with DNA [20(: the essential physics of the 
ion mediated interaction between these biomolecules is 
well captured by a simple two plates model, which opens 
the way to analytical estimates of the location and depth 
of the corresponding energy well. 

In salt-free solutions with spherical counter-ions of 
size b, the threshold distance D* MF — h* MF — b at 
which the electrostatic attraction is dominated by the 
osmotic repulsion for two plates bearing uniform sur- 
face charges o\e and 02 e (with a\U2 < 0), is simply 
given, within mean-field, by the difference of their re- 
spective Gouy-Chapman lengths: D* MF = — /x 2 | [l9j . 



These quantities read 1^1(2) = [^qIbWi(2)\] 1 , where 
Ib = e 2 / (A-TreksT) is the Bjerrum length -about 0.7 nm 
in water at room temperature- that is defined from tem- 
perature T and solvent permittivitty e. Relying on the 
Poisson-Boltzmann MF approximation, the previous re- 
sult only holds provided the Coulombic coupling between 
counter-ions is not too large. More specifically, this 
means that the two coupling parameters Si and S2 - 
defined as Hj = 2Trl B q 3 \ai\- should both be small [3. l2lj. 
However, in cases of practical interest with multi-valent 
counter-ions, the coupling parameter may be large; for 
instance, converting the charge of double-stranded DNA 
into an equivalent surface charge, one finds S ~ 23 in 
water at room temperature with di-valent ions (q = 2) 
and S ~ 76 with q = 3 (22|. In this paper, our goal 
is therefore to study the fate of the attraction/repulsion 
transition for oppositely charged interfaces, under strong 
coulombic coupling (large S limit). 

For the following discussion, it is instructive to remind 
the essential features of a single strongly coupled planar 
double-layer, without added salt (i.e. counter-ions only 
do ensure electroneutrality) |21H32| |. Irrespective of the 
value of S (from mean- field to strong coupling), the typ- 
ical distance that counter-ions may explore away from 
the charged wall is given by the Gouy length [i defined 
above. At large S values, the counter-ions form a strongly 
modulated liquid (if not a true crystal at asymptotically 
large S), with a typical distance between ions measured 
by a± = \J q/ira |23l |. as required by electroneutrality 
(ct7to^ — q). It therefore appears that /x -C a± when 
S > 1, where a± measures the size of the correlation 
hole around each ion. As a consequence, the counter- 
ions, that form a strongly correlated liquid parallel to 
the plate, effectively decouple in the direction perpen- 
dicular to the plate, and the leading order profile in the 
strong-coupling expansion is given by the interactions of 
individual counter-ions with the confining charged inter- 
face [2ll - [3"2l | : this single particle picture simply yields a 



leading exponential counter-ion density profile, with a 
characteristic length fi. Counter-ions interactions con- 
tribute to the sub-leading terms [3l| , and will not be ad- 
dressed here: we shall restrict to the low hanging fruits of 
the single particle viewpoint, that provides the dominant 
strong-coupling behaviour. We also stress that again for 
large S, we not only have /iCai but also a± <C q 2 lB- 
More precisely, it is useful to keep in mind the following 
relations 



(1) 



where the numerical constants are immaterial. 

In the following, we shall consider two uniformly 
charged plates 1 and 2, with respective charge densities 
a \ < and < 02 < Plate 1 is neutralized by 

counter-ions of valency q > while —q counter- ions neu- 
tralize plate 2. The corresponding micro- ions remain in 
the gap of width h between the plates to ensure global 
electroneutrality (see. FigJT]). Our goal is to character- 
ize the strong coupling regimes, and to infer the equa- 
tion of state at short distances from the knowledge of 
ionic density profiles, making repeated use of the con- 
tact value theorem [33l l34j , that will be reminded in due 
time. Several situations will be worked out, depending 
on the formation of Bjerrum pairs +q/ — q between op- 
positely charged micro-ions. In addition to Si and B2, 
the physics of the problem is thus ruled by another cou- 
pling parameter, T, to be introduced in section Til Al and 
that quantifies the tendency to form +qj — q pairs. This 
short range study is developed in sections IIIBI and III CI 
It will be complemented by a large distance analysis in 
section lPl Dl from which a tentative "phase diagram", al- 
lowing to discriminate repulsive from attractive regions, 
will be put forward. Conclusions will be drawn in section 
IIII and the possible relevance of our approach to weak 
couplings will be discussed. 



II. STRONG-COUPLING APPROACH FOR 
OPPOSITELY CHARGED PLATES 

A. Crowding versus pairing 

Whereas previous works pertaining to the strong- 
coupling limit have been mostly performed in the limit 
of point counter-ions, it is possible to transpose in some 
cases the results to the case of finite size ions, essentially 
by taking 6/2 (the ionic radius), as the ion-plate dis- 
tance of closest approach. For a single plate, the density 
profiles in the two cases are therefore identical, up to a 
coordinate shift z z — 6/2, where z denotes the dis- 
tance to the plate. Likewise, in the two plates problem, 
the plate-plate distance of closest approach is 6. More 
precisely, the 6 = and 6^0 cases coincide provided 
packing effects are negligible (see Fig. [JJ-left) while in- 
creasing ionic size 6 necessarily leads to a situation where 
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FIG. 1: Schematic view of the two-plates system. Micro-ions 
are hard spheres of diameter b with charges +qe or —qe. The 
width of the slab between the plates is denoted by h, and we 
define D as h — b. 




a± 




FIG. 2: Two strongly coupled double-layers in different 
crowding regimes. The left hand side plot shows the un- 
crowded situation where the finite ionic size does not perturb 
the point-counter-ion predictions, while on the other hand, 
the right hand side plot is for a case where hard core effects 
lead to crowding, with bi or multi-layers of counter-ions in 
the vicinity of the plate. For this example to be meaningful, 
we set q = +2 for the valency of the counterions. 



6 becomes of the order of a± , so that the double- layer can 
no longer accommodate a mono-layer of counter-ions (see 
Fig. [H-right). Understanding the behaviour of strongly 
coupled and crowded double-layers is an open problem 
that lies beyond the scope of the present work, so that 
we will restrict to the cases where 6 < a±, i.e. to not too 
big micro-ions. This requirement should be enforced for 
both plates: 6 < a^^K 

In addition to crowding, micro-i on p airing may take 
place in the two plates problem 0, |35[, see Figj3l The 
tendency for +q and —q micro-ions to form neutral pairs 
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Boltzmann law: 



27re(fe l + 6 2 



FIG. 3: Pair breaking/formation mechanism. Starting from 
a contact configuration for two ions of opposite charges qie 
and — q2e on the left hand side, the tendency to remain in this 
configuration, or, on the contrary, to break the pair, is given 
by comparing the electrostatic loss to the thermal energy. 



at T 7^ is quantified by the ratio between the direct elec- 
trostatic interaction at close contact and fc^T i.e. F" = 
q 2 lB/b. Interestingly, keeping in mind the "no-crowding" 
condition sketched above (b < a 1 ^'^), we get the in- 
equality 2r 2 > H (where 5 = rnax(Hi, S2) = Hi), so that 
the possible values T can take are bounded from below 
by y/a/2. Consequently, strongly coupled uncrowded 
double-layers lead to the important formation of Bjerrum 
pairs (large T). We nevertheless start by considering the 
rather narrow region where 3 J$ \/H/2 < T < 10 in which 
pair formation can be neglected. The above constraint 
translates into 20 < 5 < 200, where S is large enough to 
allow for a strong coupling analysis in due form to unveil 
the main features. A more quantitative description pre- 
sumably requires, especially at the smaller S- values in- 
volved, an intermediate approach interpolating between 
the mean-field and strong-coupling limits |36h39| . 



B. Small separation distances without pair 
formation 



We define in the subsequent analysis D as the shifted 
distance between the plates : D = h — b. The first sit- 
uation addressed is that where Bjerrum pair formation 
can be neglected, which is the assumption made in (25[. 

Under strong coupling, if D < a± (which implies that 

(2) 

D < a]_ since 02 < |ci|), the single particle picture 
where each micro-ion only interacts with both plates and 
not with its fellow micro-ions is valid. A tagged micro- 
ion feels an electric field — 27r(|<7i| + 02) .2 /e where £ is a 
unit vector along the z direction, hence a linear poten- 
tial in z. The corresponding number densities n + (z) and 
n-(z) for both +q and — q species follow then a simple 



n±(z) = n ± e 



(2) 



where n+ and are two normalization constants, and 
where we introduced the reduced distance to plate 1, 

1 1 1 



z = z/X with 



A jJLl LL 2 ' 



(3) 



The two factors 



,(()) 



can be determined from the clcc- 



troneutrality conditions: 
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so that 
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The expression for the reduced pressure 

(3P 



n = 2ttI b q 2 n\ PP 



2-kIb o{ 



(4) 
(5) 

(6) 
(7) 

(8) 



directly follows from the contact value theorem, which 
yields the pressure P in the form [33, 34]: 



(9) 



Consequently, we have 
11(D) = C coth(|)+i(l+C 2 ) 



coth f -1 



(11) 



where we introduced the charge ratio C = cr 2 / 1 o"i | . Equa- 
tion (fTTj) is independent of the plate (1 or 2) where the 
contact theorem is applied, which provides a consistency 
test for the argument. In other words, the pressure P 
is invariant under the change £ — > l/£, so that the re- 
duced pressure should change according to II — > n£~ 2 
when £ — > £ _1 . This property can be checked directly 
on Eq. More importantly, expression (jlll) is pos- 

itive for D > 0, for all values of the charge ratio, see 
Fig. [H Therefore, the interaction between two oppo- 
sitely charged plates is always repulsive at short dis- 
tances in this regime. The physical mechanism behind 
this repulsive behaviour is the following. Compared to 
the large distance limit where n + (b/2) ~ 2itIb<J 2 , as fol- 
lows from Eq. © and the fact that both the pressure and 
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n_(b/2) vanish, bringing the plates at short distances 
where D < a± enhances the electric field felt by +q 
micro-ions, which has the result to increase their den- 
sity at contact with plate 1. Invoking again the contact 
theorem ©, the consequence is that P > 0. The in- 
teractions between +q and — q micro-ions could counter- 
balance this effect, but these interactions have been dis- 
carded here, with the neglect of Bjerrum pair formation. 
We will see below that +q/—q interactions, when rele- 
vant, completely change the phenomenology. 




FIG. 4: Plot of the reduced pressure following from the short- 
distance equation of state pip , as a function of rescaled dis- 
tance D, for several values of f = <Jil\a\\: C — 0.25 (dashed 
line), C = 0.1 (dot-dashed line) and ( = 0.025 (solid line). The 
short-distance requirement D < a± (the so-called Rouzina- 
Bloomfield criterion [23]), translates into D < H 1 ^ 2 . 



C. Small separation distances with pair formation 

We now turn to the case where 1 <C < T, with a 

strong tendency for two micro-ions to form a neutral pair. 
While pair formation is unlikely as long as the two con- 
densed micro-ion layers from each plate do not overlap, it 
turns important at smaller separations. By "pair forma- 
tion" , we loosely refer here to the more or less complex 
structures, or aggregates, that may form from the asso- 
ciation of several of individual pairs. Pairs may indeed 
exist in the form of well defined entities, but may also 
self assemble into chains, see e.g. Ref. [42j |. or into more 
complex structures (regular or empty crystals) uncovered 
in a related context in Ref. [43|. The corresponding ag- 
gregates are electrically neutral, with number of Bjerrum 
pairs involved per unit area limited by the less abundant 
species of micro-ion, i.e. the counter-ions of the plate 2. 
Therefore, the aggregate surface density is bounded from 
above by <t 2 /<7. These aggregates coexist with a strongly 
correlated Wigner-like crystal made up of the remaining 
majority species. In this work, we did not attempt at a 
precise evaluation of the aggregate, or "pairs" contribu- 
tion P agg to the total interplate pressure P, but instead, 
we considered two limiting cases, where we bound P a gg 
from below by (see section III C II) , and from above by 
kTa 2 /(qD) (see section HTTT21) . The latter bound corre- 
sponds to a density of aggregates, that are neutral en- 
tities, equal to 01 /(qD), that is, to the maximum mean 
density of possible pairs. Any self-assembly of the pairs 
in a more complex architecture leads to a decrease of that 
density. We now investigate separately these two limiting 
cases. 



1. Without the osmotic contribution from the pairs 



Three remarks are in order here. 1) We see that the 
reason for observing repulsive behaviour (enhanced elec- 
tric field acting on a micro-ion within the single particle 
picture) is the same as that leading to attraction in the 
like-charged case (decreased electric field, with a corre- 
sponding decrease of micro-ionic density at contact; this 
effect is most pronounced in the symmetric case where 
<7i = for which the electric field vanishes and the 
micro- ion densities is uniform in the z direction). 2) The 
possibility of attraction under strong coupling reported 
in 25] for oppositely charged plates stems from the fact 
that only one type of micro- ion was considered in [2o| . 
This results in a smaller amount of counter-ions (per unit 
surface) , compared to that which is necessary to neutral- 
ize the isolated plate, with a concomitant decrease of 
contact ionic density, which opens the way to a possible 
attraction. The physical situation considered in [25( thus 
differs from the one under study here. 3) Our finding 
P > relics on the condition D < . At large dis- 
tances, we should have P < since a mean-field scenario 
is then expected to prevail @, Hfl H(J 53 . We will come 
back to this point in section III Dl 



With only positive counter-ions in the system, the typ- 
ical lateral distance becomes A± = y/q/Tr(\ai \ — 02) and 
for separation distances D less than A± , the single parti- 
cle picture holds, and yields the micro- ionic density n+ (z) 
in a form similar to @: 



7i+(z) 



(12) 



where n+ is a positive constant. As for n+ , the rffl 
prefactor can be determined using the electroneutrality 
condition: 



-h-b/2 

q I dz n + (z) = \a\\ — a 2 

6/2 



and the ion density then reads: 



n+ (z) 



<ri\ 



q\(l 



a2 c b/2-z 
J>-h\ 



(13) 



(14) 



In a first step, we do not consider the contribution of 
Bjerrum pairs to the total pressure (we therefore bound 
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from below the term P agg by 0). In doing so, the pressure 
at a given reduced separation distance D can again be 
found by means of the contact value theorem, with only 
one species of micro-ions: j3P = n+ (b/2) — 2-k\bo\ 1 so 
that 

11(5) = ~(1 + C 2 ) + ^(1 - C 2 ) coth(|) (15) 

We recover the same expression as in p5| . from a "me- 
chanical" (contact theorem) instead of "energy" route. 
Unlike Eq. fTTj) . Eq. (fl5)) does not have a definite sign and 
as a consequence, the interaction is attractive at large 
distances: there exists a threshold 

DS = 2Aln(H/cr 2 ) (16) 

below which the interaction becomes repulsive (II > 0), 
see Fig. [3] We note from Eq. (|16l) that attraction pre- 
vails until D — > when 02 — > \<J\ \ and we add that as long 
as micro-ions (even in small amount) remain between the 
plates as required by electroneutrality, the correspond- 
ing entropy cost for confinement makes the pressure di- 
verge (hence positivity) for D — > 0. Only for £ = 1, i.e. 
(T2 = —ci would micro- ion total density vanish, which 
leaves two oppositely charged plates interacting without 
any screening. In that specific case, the interaction is ob- 
viously attractive until close contact h = 0. On the other 
hand, in the large D limit -with nevertheless the D < A± 
requirement enforced-, and for any value of £, one can 
immediately find the pressure from the contact theorem 
applied at plate 2: for large D, the positive micro- ions 
are expelled from the vicinity of the positive plate, which 
means that the contact density n + (h — d/2) vanishes and 
that f3P — > —2-kIbo~2. In terms of rescaled pressure, we 
then have II — > — £ 2 , which is indeed observed in Fig. [5j 

2. With the osmotic contribution from the pairs 

We now include the pairs contribution to the equation 
of state, through the upper bound a-ij^qP)) alluded to 
above. We then get 

n= (1^0C_l (1 + C 2 )+ l (1 _ c 2 )coth( ^ ) 

Clearly, compared to expression (fT5)l . the effect of this 
osmotic contribution is to increase the threshold value 
where repulsion (II > 0) can be observed. The two lim- 
iting behaviours, Eq. (fT5")l and Eq. (|T7|) . are sketched in 
Fig. [5] The corresponding values of the thresholds Dq 
and D\ are indicated. These two quantities are plotted 
in Fig. [S] as a function of charge asymmetry, together 
with the analytical estimation of D\ obtained as follows. 
If D is large enough, Eq. (|T7| simplifies to: 

n--e 2 + (1 t c)c . (is) 
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FIG. 5: Pressure curves from Eqs. (|15[) (solid line) and (|17|) 
(dashed line) for Q = 0.2. 



This expression can now exhibit a repulsive behaviour 
below 5 = (1 + 0/C i-e. 

D{ ~ M 2 (19) 

It can be seen that this approximation (dashed line) is in 
fair agreement with the root of Eq. (jTTJ) found numeri- 
cally (dotted line), in the whole available range. 

The single particle picture invoked here relies on two 
assumptions. First, T should be large, to have pair forma- 
tion. Second, the (shifted) distance D between the plates 
should not exceed A±, see the upper dot-dashed line in 
Fig. [5] Making use of approximation (TT91 . this means 
1/(2S2 + 1) < C- F° r smaller values of £, the analysis 
is significantly more complex (loss of the single particle 
view-point). Hence, if ( > 1/(2^2 + 1) the interaction is 
repulsive at short distances and then turns attractive at 
intermediate distances, while if ( < 1/(2S 2 -I- 1), Eq. (fTS|) 
does not lead to any transition between repulsion and 
attraction, and is always repulsive in its range of valid- 
ity, as was the case for Ea. lfTTj) . Of course, for large 
S 2 , the threshold (2S2 + 1) _1 is small, so that extremely 
asymmetric cases only (very low £ are not covered by our 
analysis. 



D. Large separation distances 

Our analysis has so far been restricted to short dis- 
tance expansions. We are now interested in large dis- 
tance asymptotics and in attempting to match the short 
and large distance behaviours. In doing so, we will dis- 
cuss qualitatively an attraction/repulsion transition of an 
"effective" mean-field type, which leads to reentrant at- 
traction as the distance between the two plates is varied 
from infinity to close contact. 
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FIG. 6: Reduced threshold distance discriminating (strong- 
coupling) attraction from repulsion, as a function of £, in 
semi- log scale. The solid line shows Dq obtained from (|15[1 . 
the dotted line is for D1, the exact root of Eq. (|17[) . and the 
dashed line displays the approximation (|19|l . The upper dot- 
dashed line shows Aj_, in the particular case E 2 = 25: our 
approach is meaningful for D < A± only, which translates 
into C > 1/(2E 2 + 1), see text. Here, 1/(2H 2 + 1) ~ 0.02. 



1. Crossover between strong- coupling and mean-field 
regions for one plate 

We will assume first that a given strongly-coupled plate 
(having thus a large Sj), can be effectively described by 
mean- field theory, at sufficiently large distances z. This 
"common wisdom" stems on the remark that for large 2, 
the typical distance between counter-ions becomes large, 
which leads to a low coupling regime [H, [H, H(J l4lj . 
It should be emphasized though that the above point of 
view, that predicts a large z density decay in l/z 2 , is 
incorrect in two dimensions, as shown in a recent work 



31]. The present study pertains to three dimensional 
systems, so that we nevertheless expect for a single plate 
the crossover scenario discussed in Ref . [4lj ] , and summa- 
rized in Fig. [7] In essence, the density is expected to de- 
crease exponentially at short distances, and algebraically 
at large distances: beyond a distance 5 from the plate, 
the counter-ion density n is simply given by the solution 
of the nonlinear Poisson-Boltzmann (PB) equation: 



1 



2vr/ s q 2 (3 + ^ ff ) 2 



(20) 



where 3 = z — 6/2 and where li is an effective 
Gouy-Chapman len gth characterizing this long range be- 
haviour. Following [4ll ] . one can match the two regimes 
by assuming that the condensed counter-ion layer forms 
a 2D One Component Plasma and by applying a mean- 
field approximation for the dilute layer. Equating the 



two corresponding chemical potentials yields 
n(S) =n sc e~^ 



(21) 



where 

TIqq is cl ^-related average density in the condensed 
layer and /3e c (2) w — 1.56\/S is the contribution to the 
2D one component plasma chemical potential that stems 
from the correlations between the counter- ions @, |24 |. 
Extrapolating the validity of Ea. (|2"0"|) to 5 — > 6/2 and as- 
suming that in such a situation, n sc is well approximated 
by the avera ge d ensity over the characteristic length l sc , 
we arrive at 1411: 



" (s = 0) = s; 



1 



-l at ./fl 



-1.56VH 



(22) 



Equation (|2"2"j) is nothing but the density that the mean- 
field profile, valid at large distances from the plate, would 
have if extrapolated at 3 — 0, and is therefore not the real 
density at the plate. However, invoking Eg . (|2"0")) , it al- 
lows one to estimate the effective Gouy-Chapman length 
ll cS corresponding to the charged plate dressed by a con- 
densed counter-ion layer, which will prove useful in the 
following. For the subsequent quantitative discussion, we 
shall take the value l sc ~ 3.6/x, already used in [4lj . 



Counterion 
Density 




PB 



FIG. 7: Crossover scheme: following Ref. [4l]], schematic rep- 
resentation in semi-log scale of the crossover between strong 
coupling (SC) and mean-field regimes for the counter-ion den- 
sity in the vicinity of a single highly charged plate. The solid 
curve represents the exponential decay expected close to the 
plate while the dashed curve stands for the algebraic decay 
expected at large distances where the Poisson-Boltzmann the- 
ory (PB) should hold. The shaded region corresponds to the 
crossover between these two regimes: l sc is the distance be- 
yond which the exponential decay no longer holds and S is the 
distance to the plate beyond which the PB profile is expected 
to be valid. 



2. Application to the two plates problem 

When the separation distance between plates 1 and 2 
is decreased from infinity, the first inter-plate weak in- 
teraction regime is expected to be of mean-field type, so 
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that the presumably large distance attraction may turn 
into repulsion at a distance D* MF = Ifif 1 — In this 

picture, the distance is varied at constant effective Gouy- 
Chapman lengths /if 1 and /i^ given by: 

fif = 1.92 /i ie °- 78VH *, i = 1,2. (23) 

If |/i° ff — /x| | is significantly larger than the char- 
acteristic thresholds obtained in the previous subsec- 
tions, we should have the following "reentrant" sequence 
{attraction — > repulsion — > attraction — > repulsion} as D 
decreases. The first transition is described by a mean- 
field argument, and the last one by strong-coupling con- 
siderations, but the intermediate transition {repulsion 
— > attraction} occurs in a crossover region that re- 
sists our theoretical understanding, and where additional 
{repulsion — > attraction} transitions might take place. A 
related question deals with the lower bound for the dis- 
tance Z?bound below which the mean-field profiles are no 
longer accurate. For the sake of completeness, we will 
consider below that -Dbound = Depending on 

the respective surface charge densities a\ and 02 , we can 
then discriminate between two distinct situations: 

• |^° ff — < D bound- The interaction between the 
two plates is always attractive at large distances 
(mean-field regime), and then at short separation 
distances, the strong-coupling phenomenology de- 
scribed in the first part of the paper prevails. 

• I Mi — Ml I > Di, oun d- There is then already a 
transition between attraction and repulsion in the 
mean-field regime. By decreasing further the dis- 
tance D and entering the short distance limit, one 
should observe another attractive range, as ex- 
pressed in Ea. (fT5|) for instance, before repulsion 
sets in at even smaller separations. 

More complicated scenarios could be envisioned, but we 
summarize in Fig. [5] the simplest possible, and provide a 
phase diagram obtained when considering that the Bjcr- 
rum pairs do contribute to the pressure, as in section 
III C 21 We note that, for the parameters chosen, there is 
a reentrant behaviour observed with respect to the sep- 
aration distance, in a large fraction of the (£, D) plane, 
more specifically when £ < 0.5 (this threshold depends 
on the value of S 2 chosen, and increases with S 2 ). We 
remind that the bottom part of the diagram, more specif- 
ically for £ < 1/(232 + l)j corresponds to a region where 
our arguments do not apply, as discussed in section lll C 21 
In this region, our short scale analysis provides an "all 
repulsive" behaviour, and we may then speculate that 
repulsion persists up to the effective mean-field thresh- 
old indicated by the asterisks, which corresponds to large 
distances, on the order of 100<7 2 ^b or more. 

III. DISCUSSION AND CONCLUSION 

In this paper, we have analyzed the interaction of two 
oppositely charged parallel interfaces, each neutralized 
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FIG. 8: Sketch of attractive and repulsive regimes, as a func- 
tion of the ratio £ = (T2/\(Ji\ and the distance D (in log 
scale). The repulsive island on the left hand side, delim- 
itated by a continuous line, shows D\ where the pressure 
in Eq. (|17[) vanishes. The other repulsive region on the 
right hand side, delimitated by a dotted line (star symbols), 
shows D* MF = |/xf fl ' — /i2 ff I, where the effective Gouy-Chapman 
lengths are given by Eqs. H23[) . These data are displayed pro- 
vided they satisfy the constraint D* MF > Dbound = ffl'J' +0^. 
Likewise, the left boundary for this mean-field repulsive island 
has been taken to be Dbound, that is shown with a dotted line 
(square symbols). Here, we have taken H2 = 10. 



by its own counter-ions, without other micro-ions in- 
volved (salt-free case, but with two species of micro- 
ions, of opposite signs). We have shown that a repul- 
sive behaviour, quite expectedly, is always present at 
short enough separations; it simply stems from the di- 
verging entropy cost for confining micro-ions in a slab of 
vanishing extension. Our analysis completes the known 
Poisson-Boltzmann phenomenology by investigating the 
case of strong coulombic couplings. Short distance ex- 
pansions reveal that depending on the formation of Bjcr- 
rum pairs between the oppositely charged micro-ions, an 
attractive regime may or may not be observed. By for- 
mation of pairs, we understand here the wealth of differ- 
ent self-assembly scenarios where the pairs may further 
associate into more complex objects, such as chains or 
various crystals [13, El] . We did not attempt at a pre- 
cise evaluation of the corresponding contribution to the 
pressure -a particularly demanding task- but rather, we 
analyzed limiting cases where this unknown contribution 
is bounded by reasonable values, see section III CI We 
have supplemented our short distance analysis by a more 
speculative investigation of the large distance behaviour, 
from which a "phase diagram" was put forward, with 
reentrant features between attraction and repulsion as 
the distance D between the plates is varied. The exper- 
imental observation of such features would imply that 
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other sorts of interactions, such as van der Waals, do not 
modify the main effects uncovered. 

In our approach, single particle arguments play a cru- 
cial role and allow us to compute the density of micro- 
ions, from which the equation of state follows. These 
single particle arguments, however, are a priori not re- 
stricted to strongly coupled interfaces, but can equally 
be invoked when the coupling parameters Si and S2 are 
small (see e.g. section 3.3 of Ref 01 an d in particu- 
lar Fig. 16 for simulation results backing up this state- 
ment in the like-charge case o~i — 02 )• Indeed, when 
D becomes smaller than the characteristic lateral dis- 
tance aj_ between ions, these ions effectively decouple 
and feel the external potential of the plates only (we 
are here concerned with the ionic density dependence 
on the z coordinate, perpendicular to the plate; in the 
transverse direction, parallel to the plate, a correlation 
hole remains around each particle, of typical size q 2 ls)- 
As a consequence, the pressures given by Eqs. dill) . 
(fT5j) and (|17l) still hold with the same range of valid- 
ity, for Si — > 0. In the corresponding distance range, the 
Poisson-Boltzmann results break down due to discrete- 
ness effects [it is therefore essential here to make a clear 
distinction between Poisson-Boltzmann theory, and the 
low S limit of the original model dealing with discrete 
particles ; Poisson-Boltzmann approach considers from 
the outset continuous density fields, and can therefore 
not be expected to hold at separation distances such that 
discreteness effects do matter, i.e. when D < a±; the ade- 



quacy of Poisson-Boltzmann to describe the low S physics 
should then be understood as a statement which excludes 
a small range of short separations D]. From the analysis 
of section HH we learn that when Bjerrum pair forma- 
tion can be neglected, the threshold distance D* where 
repulsive behaviour sets in is still given by the Poisson- 
Boltzmann result |/X2 — Mil = A*2 — Ml) provided this length 
is larger than both characteristic distances and a±\ 
In the opposite case, when |/i2 — Mil < in£(a± , a^), we 
may speculate that D* lies between aS^ and a± , since 
the single particle argument which holds at smaller sep- 
arations leads to repulsion, while the Poisson-Boltzmann 
theory yields attraction at larger separations [20] . If on 
the other hand, Bjerrum pairs form and contribute to the 
pressure through their mean density (see section Hi C 2[) . 
we have seen that D* = /i2, which is thus larger than 
the Poisson-Boltzmann result M2 — Mi- However, this re- 
sult only holds provided £ > (1 + 2S 2 )~ 1 ~ 1 (we are 
still considering the low 3j limit). Given that £ < 1 
by definition (i.e. 02 < |ci|), we see that here, the sin- 
gle particle picture does not apply up to y,i (except in 
a small £ region close to 1), which means that D* is 
larger than Aj_ = [7r(|(Ti| — o^)] -1 / 2 . It can be checked 
that generically, this length is smaller than the Poisson- 
Boltzmann prediction /i2 — Mil except again in a small 
(^-region around unity. 
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